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A string-inspired effective theory of gravity, containing Gauss-Bonnet invariant interacting with 
a scalar field, is considered in view of obtaining cosmological dark energy solutions. A Lagrange 
multiplier is inserted into the action in order to achieve the cosmological reconstruction by selecting 
suitable forms of couplings and potentials. Several cosmological exact solutions (including dark 
energy of quintessence, phantom or Little Rip type) are derived in presence and in absence of the 
Lagrange multiplier showing the difference in the two dynamical approaches. In the models that we 
consider, the Lagrange multiplier behaves as a sort of dust fluid that realizes the transitions between 
matter dominated and dark energy epochs. The relation between Lagrange multipliers and Noether 
symmetries is discussed. 
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I. INTRODUCTION 

Astrophysical data indicate that the observed universe 
is in an accelerated phase This acceleration is in- 
duced by the so-called dark energy (see Ref. [2| for a 
recent review and references therein) which nature and 
properties are not yet understood at fundamental level. 
In the most theoretical models considered in the litera- 
ture, the dark energy is constituted by some ideal fluid 
with a specific equation of state (EoS) sometimes exhibit- 
ing non-standard properties like negative pressure and/or 
a negative entropy. On the other hand, dark energy can 
be considered as a global phenomenon associated with 
modifications of gravity Q. In fact, its presence could 
point out that Einstein's General Relativity cannot be 
retained as the final theory of gravitational interaction 
at cosmological scales. In this sense, General Relativ- 
ity presents problems at UV (Quantum Gravity) and IR 
(cosmology) scales [4J. A similar situation appear for 
dark matter phenomena. No fundamental candidate has 
been revealed up to now and dynamics of self-gravitating 
structures could be addressed by modifications of gravity 
1- 

A further problem is that it is not clear why dark en- 
ergy had no effect at early epochs while it gives domi- 
nant contribution in today observed universe. According 
to the latest observational data, dark energy currently 
accounts for about 73% of the total mass-energy amount 
of the universe (see, for example, Ref. @). 

The main feature of dark energy is that its EoS pa- 
rameter wd is negative: 

wd = Pd/pd < , (1) 

where pd is the dark energy density and p^> the pressure. 
According to the standard cosmological model, this prop- 
erty gives rise to the reported apparent acceleration of the 



Hubble fluid. Although current data favor the standard 
ACDM cosmology, the uncertainties in the determination 
of the EoS dark energy parameter w are still too large, 
namely w — — 1.04j^Q < j ) o- Hence, one is not able to deter- 
mine, without doubt, which of the three cases: w < — 1, 
w = —1, or ro > —1 is the one actually realized in our 
universe @, Future observations should better con- 
strain this range of values giving also indications on the 
nature of dark energy Q . 

In order to explain dark energy at very fundamen- 
tal level, string/M-theory could suggest reliable effective 
models to be compared with observations. In particular, 
further gravitational terms, emerging from this theory, 
could become important at current, low-curvature uni- 
verse (being not essential at intermediate epochs from 
strong to low curvature). For instance, in the study 
of string-inspired gravity near the initial singularity, the 
role of the Gauss-Bonnet (GB) topological term, coupled 
with scalar potentials, is important for the occurrence 
of non-singular cosmology [ToL ITU). The dilaton coupled 
to higher-order curvature corrections, for example, as- 
sumes an important role near the initial singularity, as 
discussed in [T^ |. Furthermore, string-inspired gravity 
with Gauss-Bonnet term interacting with a scalar field 
has been proposed as a realistic candidate to address the 
dark energy issue [T^|. Specifically, as it was shown in 
Ref. fl3l| . Gauss-Bonnet dark energy may lead to the oc- 
currence of phantom cosmology without ghosts: in this 
case, the dilaton is a canonical scalar. Further aspects 
of Gauss-Bonnet accelerating cosmology have been dis- 
cussed in detail in [T3, EH- 

Recently, a new dark energy model has been proposed 
fl6L . It consists in considering two scalars where one 
of them is given by a Lagrange multiplier. This multiplier 
puts a natural constraint on the form on the second scalar 
field (in particular on its coupling and self-interacting po- 
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tential) allowing that the emerging dark energy behavior 
evolves in a dust-matter dominated era, as requested by 
observations going back in the redshift z. The interest- 
ing feature of this approach is that the whole system 
contains a single dynamical degree of freedom and this 
fact allows to solve several shortcomings related to the 
fine-tuning of ACDM model, among them the cosmolog- 
ical constant problem [23]. It is important to stress that 
such a Lagrange multiplier technique can be related to 
the existence of Noether symmetries and then it is a gen- 
eral approach to reduce dynamical systems and find out 
exact solutions, as we will discuss below fl8l . 



II. SCALAR GAUSS-BONNET GRAVITY WITH 
LAGRANGE MULTIPLIER 

Let us study accelerating cosmology in string-inspired 
scalar Gauss-Bonnet gravity with Lagrange multiplier. 
To this end, a suitable action has the following form: 
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e{!p)G - A ( -d^ifi + E%) 



(2) 



The extension of f(R) gravity via the addition of 
a Lagrange multiplier constraint has been proposed in 
Ref . |19| . Such model can be considered as a new ver- 
sion of modified gravity because dynamics and cosmo- 
logical solutions are different from the standard version 
of f(R) gravity without such constraint. This result is 
clear from a dynamical viewpoint: Lagrange multipliers 
are anholonomic constraints capable of reducing dynam- 
ics (20l. l2lj|. Furthermore, using the Lagrange multiplier 
approach helps in the formulation of covariant renormal- 
izable gravity |2^ j. 

In the present paper, we study the Gauss-Bonnet grav- 
ity with Lagrange multiplier constraints in view to re- 
cover realistic dark energy behaviors. Technically, we are 
considering a scalar-tensor-higher-order gravity where a 
Lagrange multiplier is considered. We explicitly show 
that it is possible to derive new cosmological solutions 
in this context. In particular, we compare the acceler- 
ating solutions in string-inspired Gauss-Bonnet gravity 
with and without the Lagrange multiplier term. It is 
demonstrated that, in the version with Lagrange mul- 
tiplier, one gets large number of new accelerating cos- 
mologies, including the phantom cosmologies where the 
dilaton kinetic term is canonical. As a more interest- 
ing example, the so-called Little Rip cosmology can be 
recovered. 



The plan of the paper is the following. In Sect. II, 
we discuss the scalar Gauss-Bonnet gravity in presence 
of a Lagrange multiplier. In particular, we consider the 
Friedmann-Robertson- Walker (FRW) cosmology and de- 
rive exact solutions according to the choice of the scalar 
field ip. The general scheme of cosmological reconstruc- 
tion of scalar Gauss-Bonnet gravity with Lagrange multi- 
plier is pursued in Sect. III. Here we work out the whole 
method starting from the cosmological equations up to 
the final cosmological solutions. In Sec. IV, we rewrite 
the cosmological Friedmann equations as an autonomous 
system of first order differential equations and study its 
critical points. Sect. V is devoted to the summary and 
discussion of the results. The general discussion of the 
Lagrange multiplier method in view of the Noether Sym- 
metry Approach is reported in Appendix A. 



Here A is the Lagrange multiplier scalar, e, uj and U 
are arbitrary scalar functions and G is Gauss-Bonnet in- 
variant is 



AR^BT+FP. (3) 



The above effective action represents a string-inspired 
gravity which has been mainly considered for exponen- 
tial potentials and without last term (for a single scalar) . 
The interpretation of above action is related to specific 
compactification of superstring theory in four dimensions 
where, apart from the dilaton scalar field, one more scalar 
(given by the Lagrange multiplier) is considered. 

The gravitational field equations are derived by varying 
with respect to the metric and assume the following form: 



2k 2 



+V{<p)+\[-d p <pd>'<p + U(<p) ) \ + 



(4) 



+4( _ R a0 g^)VcNpe{ip) - m a l y a v v e{i P ) + 

+2EV (U V„e(^) . 

Let us now consider a FRW metric with flat spatial 
part: 



ds 2 = -dt 2 + a{t) 2 ( dxi 



(5) 



i=l,2,3 



Then by the variation over A, we obtain the further equa- 
tion 



The cosmological Friedmann equations are then 
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de 

(oj(ip) + 2\)U(p) + V(<p)+24H 3 -, (7) 
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dt 



(8) 



at 



It is straightforward to derive, from the above equa- 
tions, the effective EoS-parameter w as: 



u}(ip)U{tp) - F(<£>) - 8iJ 



2 (T£ 
dt 2 " 



(u( l p)+2\)U( l p) + V( l p) + 24H3f t 

The equation obtained by varying the action over the 
scalar field <p (the Klein-Gordon equation) is a conse- 
quence of second Friedmann Eq.([5]) and then is easily 
achieved. 

Eqs.([7][H]) can be easily solved. Indeed, there are un- 
known functions: H, V, if, A, u> and e which may be 
used to satisfy two equations. Therefore, if we want to 
obtain a closed system of equations, we must fix some 
of the functions V, A, to and e or select the field ip and 
metric a, as well as two of the functions V, A, uj and e. 
For example, one can select: 



w(vO = l, [%) = —. 



(9) 



Here m is a constant with the dimension of mass. Also, 
the canonical scalar is considered. Then Eq. ^ gives 



m 2 t . 



(10) 



In this case, one can express the potential V and the 
Lagrange multiplier A in terms of the Hubble rate H and 
the function s: 



X 



1 



H' 



8n 2 H i(2H - H 2 ) + He - 2H 



H z \ — - 8e 



(11) 



On the other hand, one can select the scalar field poten- 
tial and the function e, following string-inspired consid- 
erations: 



V = A + V e 



blip 



e = e e b2tp . 



(12) 



Here b\ and &2 are constants. By analogy, the Lagrange 
multiplier is chosen in the same form 



A = Ai + Aoe' 



'3V 



(13) 



with 63 being some constant. 

Let us consider now some explicit cases to achieve so- 
lutions for the above dynamical system. It is clear that 
this is a sort of "inverse scattering procedure", that is 
by fixing the form of the scalar field, we solve dynamics. 
Usually, it is exactly the contrary: by fixing scalar field 
potentials and couplings, field equations gives the scalar 
field form. 



A. The case ip ~ t 

This choice of scalar <p corresponds to a constant U(ip). 
In addition, we assume that uj = ±1. Hence, 



ip = cp t, w(</j)=7 = ±1, 



(14) 



where ipo is some constant. 

One can solve the first Friedmann equation for the 
Hubble rate. In this case, the second equation becomes 
constant. We get: 



H = - 



-ipabit 



1 



24( / 3 eo^2'« 2 y(l-2a+ y/4a(a - I)) 1 / 3 
+ (l-2a+ v /4o : (a-l)) 1/3 + l), (15) 

here a = 72<^£^K 6 e 2 *' oM (^(7 + 2A) + 2V) . 

It is easy to see that a must be either less than zero 
or greater than one. That is ^0(7 + 2A) + 2V < or 
<^2( 7 + 2A) + 2V > 1. 

Choosing the constant case, one can easily meet these 
requirements. If we assume that a is much bigger than 
one, then 



H 



-b 2 (p 



l/e. 



(16) 



For this case, it is easy to solve the Friedmann equation. 
Let the Hubble parameter is given by 



H = hoe 



ipoht 



(17) 



This is the so-called Little Rip cosmology (25| which 
can be realized also in the modified gravity context (see 
Refs.[13|). The effective EoS parameter of this model has 
the phantom nature, that means 



' J eff 



= -1 - 



2<pah 
3h 



-tpoht 



(18) 



In the Little Rip cosmology, w approaches —1 suffi- 
ciently rapidly, so that it is possible to have a model in 
which pde increases with time, but in which there is no 
finite-time future singularity. The characteristic feature 
of such a cosmology is disintegration of bound objects at 
finite time (like in the case of Big Rip cosmology) . 

As the universe expands, the relative acceleration be- 
tween two points separated by a comoving distance I is 
given by Za/a, where a is the scale factor. An observer 
sitting at comoving distance I away from a mass m will 
measure an inertial force on the mass of the order 



(19) 



^Inor = mid I a = ml lH-\-H 



Let us assume the two point force is bounded by the 
bounding force F . If F iner is positive and greater than 
fo, the two particles become unbound. This effect is the 
"Rip" produced by the accelerating expansion. 
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Eq. (fiH]) shows that Rip occurs when either H diverges 
or H diverges (assuming H > 0) [25|,[26j]. As we see, it is 
possible for H and for F incr , to increase without bound 
and yet not produce a future singularity at a finite time. 
This is the essence of Little Rip. 

From Eq. (pU for (fTfjl. we have 



F inor - ml (/wo he^ ht + h 2 e 2 ^ ht ) 



(20) 



which is positive and unbounded. Thus, -Finer becomes 
arbitrary large with increasing i, that results in a Little 
Rip. The parameter Lp h can be estimated (25[. Let us 
assume 



ip h 



where /3 ~ 3.46 x lO~ A Gyr 



(21) 



In addition, the present value of the Hubble constant is 
Hq 1 = 13.6 Gyr^ 1 . In this case, one can estimate the 
time required for the disintegration of an object of the 
size of the Solar System as t ~ 7750 Gyr. Furthermore, 
it is easy to get the theory constants values (poh = 2.99 x 
10~ 3 , ho — 70.59 x 10 -3 . For such a set of constants, the 
EoS parameter is equal to —1.0272 and this value does 
not contradict the observational bounds. 

Let us make some general comments about the impact 
of the Lagrange multipliers on Fi ner . From (fl~9l) and ([7|) 
one finds 



Indeed, we have 

H = h e voht , 
1 



V = 2^o7- 

£ = 0, 

A = - 7 - 



2ip hh ht 



^fe 2 ^ 4 , (24) 



2h Q h 

2* 



,<poht 



Once again, the Friedmann equations are satisfied. Ob- 
viously, this scheme may be applied to generate new cos- 
mological solutions. It is clear that the potential and the 
Lagrange multiplier have now the form 



fll 



X 



(25) 



This solution has been obtained in ref.(l3| in the same 
way. 

As an example one can keep just a Lagrange multiplier 
(V = and e = 0) in the form (p~3]l . Then solving the 
Friedmann equations, it is easy to find the Hubble rate 



H 



V6 



tan 



f0\/l K (t — const) 



(26) 



H z + H 



(22) 



1 



1 



--7C/ --XU+-V + AH e + 8HeH + AH i 



It can be seen that, depending on the sign, A can either 
enhance or decrease F iner . In other words, for realistic 
models, the role of Lagrange multiplier may be related 
to the increase/decrease of the time left before the disin- 
tegration of bound objects in the universe. 

Let us now find the solution of the Friedmann equa- 
tions (05]) for the model ([T71) asuming the exponential 
scalar functions (|12j|13j) . There are two solutions: 



Vi 
£1 
Ai 



1 2 H 2 
1 



or 



v 1 2 



H 2 



1 



Atpgh hoK 2 H 
2 

-7 - 



(fl 2 K 2 



H z 



or 



or 



£2 = 
A 2 = 



1 1 

16k 2 H 2 
3h 



-7' 



, (23) 
H. 



In the first case, a tends to a constant value, which 
corresponds to the solution (|T5l . In the second case a — > 
0. 

Now we consider the case when e = 0. This solution 
is studied in Ref. jl5j . Selecting the potential and La- 
grange multiplier as (p~2|) and (fl"3f . we do not get the so- 
lution in the form H ~ e v . However, it is easy to choose 
the potential in a slightly different form to get again the 
exponential Hubble rate. 



Since the solution is derived from the second Fried- 
mann Eq. ([5]). there is no dependence on the Lagrange 
multiplier. Now using a Lagrangian multiplier, one can 
convert the first Friedmann Eq. ([7]) into an identity. In 
this case, we obtain A as follows 




tan 



W ^Vo\/7' c (*- const) 



(27) 



For this solution, EoS parameter w lies in the range of 
— 1/3 to infinity. In other words, the expansion is decel- 
erating. 

There is another solution, if 7 = —1 (in this case, 
H =const). The presence of the Lagrange multiplier 
leads to new solutions which do not exist without it. For 
example, using the conditions (fl"2j) . ([M)) and (fl~7j) one 
cannot get a solution without a Lagrangian multiplier. If 
Lagrange multiplier (|13p is taken into account then it is 
easy to construct the solution (I21H) and This means 
that the presence of the Lagrange multiplier effectively 
changes the dynamical system giving rise to a constrained 
dynamics with, in principle, different solutions. 



B. The case ip ~ \nt 

This form of scalar field allows us to choose the Hubble 
rate as \/t for the functions V, e and A given by Eqs. (fT2"|) 
and ([13]). 
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That is 

H=^-, <p = Lp ]n^, 7 = ±1, (28) 

where ho, tpo and ii are some constants. When ho > 0, 
we have a quintessential power-law expansion. For 

H=~^-. tp = <p ln t -^, u(tp) =7 = ±1, (29) 



we have a phantom power-law expansion when /ig < 0. A 
phantom model describes the Big Rip finite-time future 
singularity: for t — t s the scale factor tends to infinity. 
Now we can write the Friedmann equations as 



A 



Vo7 . ip 2 \ 



2t 2 k 2 t 2 t 2 
^ 2 a x t- 2+b n^ - t b H, bl V = 

i ^7 2/iq , 3h 2 t 



- - 24e b 2 h 3 t- 4+b Hi b2 - 



(30) 



2t 2 k 2 t 2 k 2 t 2 
+ 8e b 2 h 2 (3 - b 2 - 2h ) r^t^ 2 - t b H^ bl V = 0. 

For ho > (for ho < we have to replace t with t s — t) 
and 

V = A + V Q e blip/vo , e = e e b2v/ipo , 

A = Ai + X e b3v/v ° (31) 

First, we consider the case when the Lagrangian mul- 
tiplier and the potential is zero (A = V = 0). It is easy 
to see, that 



, 6fag(l - ftp) 
fe 2 =2, ip a = -— 



7 k 2 (1 - 5/i ) 



(32) 



(3/tp - 1) 
8h {5h — \)n- 



(ts - t) , h < 



{3h - 1) u ,2 



8h Q (5h - 1)k 2 



(*) , > 0. 



The solutions exist if 

a) 7 <0 and h <0 or 0</i <l/5 (or h >l), 

b) 7 >0 and l/5<h <1. 

In the case in which the scalar field ip is canonical 
(7=1), and there is no potential (V(<^)=0), we cannot 
obtain the effective phantom cosmological solution with 
w < — 1. A similar situation is obtained for e = and 
V 7^ 0. In this case, the potential has the form 

h (3h a - 1) 
k = — , for ft > 



or V = 



kH 2 
hofihp - 1) 

K 2 (t s -t) 2 



, for fto < 0. (33) 



In addition there are restrictions on h : tpg'y = 2/i /k 2 . 
One sees that w < — 1 only in the case 7 < 0. 



Now let A = but V 7^ and e 7^ 0. In this case, the 
solution exists only for 6 2 = —b\ = 2 [l3j]. Then 



V 



(2fe - <p 2 ol k 2 ) t 2 
16h 2 (l + h )k 2 ' 

6(-l + /i )fe§ + ^7(-l + 5/i )fc 2 

2(1 + h Q )k 2 t 2 
ho > or t — >• t s - i, if /i < 0. 



(34) 



It is evident that there is a phantom solution for 7> 0. 

If A 7^ and V — 0, we obtain the following conditions 
restricting /iq 



2h 2 (3h 



(35) 



Therefore, to get the value of w about —1, 7 must be 
negative. If 7 > 0, then ho has to lie in the following 
range < ho < 2/3. For this case, we have 



3-2hn 



A = 



7(1 - Sh ) \2eoih 2 {2ho - 3) 



3h - 2 (3h - 

for ho > and t —$-t s 



l-2/io 



- t if h < 



(36) 



For sufficiently large ho, the function e has a negative 
power and quickly decreases with time (for/ip > 0) or 
increases when approaching to t s (for /io < 0). For the 
model without a Lagrangian multiplier (|33p . function e 
grows for /iq > and decreases for ho < 0. 
If A 7^ and e = 0, we get 



A = —7 - 



2h v = ~4fe + 6h 2 + ^ 2 7 K 2 (3?) 



lfi 2 K 2 



2kH 2 



Such behavior is similar to the one in the model without 
the Lagrange multiplier (f3"3"|) . 

Suppose now that the action has the form @, that is 
V 7^ 0, e 7^ and A 7^ 0. From Eqs.([30D, one sees that 
there are several solutions: 

1) b\ = — 4 + 62- For this case, there are two additional 
conditions: 62 = 2 and b 2 =/= 2. 

In the first case it is easy to find: 



V = 



32e (l - 2h Q )h 2 K 2 + (2ho(3h - 2) + <pljK 2 )tj 



2n 2 t 2 t 2 



t 2 



(38) 



. -16e ^(l + h )n 2 + (2h - ^h^t 2 _ 

A — 7y — - — 7^ — const. 

tp Z 0K 2 t{ 

This means that the solution is similar to the solution 
for the case without Lagrange multiplier (1351) . The time 
dependence is the same. The difference enters via an 
arbitrary factor So by choosing the Lagrange multiplier. 
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Now, let b 2 7^ 2, then we obtain the following condition 



± 



V2 h (2-3h ) 



7 



(39) 



Again limits appear on the Hubble rate. Phantom cos- 
mology is possible only for 7 < 0. If 7 > 0, then ho must 
be greater than zero and less than 2/3. For this case, we 
have 



V = — 



8e b 2 ho(b2-S + 2ho)(t/t 1 ) 



-4+f) 2 



£0- 



7 



1 - 3/tp 4g 6 2 7(b2 



h )h K 2 



(2 - 3/i )*i 



(40) 



-2+62 



ho > 0) or another case where the potential decreases as 
1/t 2 , and e increases as 2/l o (f or ^ < q) or decreases 
(for /io < 0). Let us consider the case with/without the 
Lagrange multiplier but with the same scalar potentials. 
Let 61 = -2, b 2 = 2 and ho = -80/3, K = 1, 7 = 1. In 
this case, we numerically derive 

Vo = 2304.8 + 4.23377 A 
e = 0.000189478 + 6.84862* 10~ 6 A. (43) 



The Big Rip cosmology occurs with the corresponding 
Rip time for model (|2"9")l at t s = 376.27 Gyr. In summary, 
the presence of Lagrange multiplier may lead to the gen- 
eration of new cosmological solutions if compared with 
scalar Gauss-Bonnet gravity without such term. The 
general reason for this different behavior will be discussed 
below. 



2) Assume now the following conditions b 2 — 3 — 2h . 
Then 



-4/i + 6h 2 + <p 2 -yn 2 



V = 



£ = E 



A = —7 - 



2nH 2 
3-2/10 



2/ln 



v?6 



2Z~2 



(41) 



2A So hl(-3 + 2h ) (£) 



1— 2hn 



III. THE RECONSTRUCTION OF 
SCALAR-GAUSS-BONNET GRAVITY WITH 
LAGRANGE MULTIPLIER 



In this section we discuss the general scheme of cosmo- 
logical reconstruction in the scalar-Gauss-Bonnet gravity 
with the Lagrange multiplier. Let u)(ip) be 1. Friedmann 
Eqs. ([7][5]) can be written as 



3) There is one more case: b 2 — 4. Then 



V 



e = 



32e hl(l + 2h ) -4h + 6h% + 



A = -7 + 



+ 



2kH 



2-/2 



2h 32e /ig(l - h )t 2 



(42) 



\h 2 = 

= (1 + 2A ( ^))^ + W + 24^«11 
2 K 2H + 3H 2 )=^-V(^-8H 2d2£iip{t)) 



(44) 



1 



dt 2 



dt 



dt 



(45) 



Hence, the Lagrange multiplier presence helps to gen- 
erate new cosmological solutions. For models without 
such multipliers, we find that the potential is propor- 
tional to the square of the time, while the scalar function 
of the Gauss-Bonnet invariant is inversely proportional 
to the time. In the model with the Lagrange multi- 
plier we get new solutions ([3T|). (gO}, (0TJ) and |g2J. In 
all these cases, the dependence of the coefficient of the 
Gauss-Bonnet invariant is explicit and simple, whereas 
the scalar potential often shows the same time depen- 
dence (V ^ 1/t 2 ). It should also be noted that there 
may occur phantom cosmologies for the canonical scalar 
field ip. Only for solutions (|4T))) we obtain restrictions on 
ho (for V, £ and A /0). 

Without the Lagrange multiplier, the scalar potential 
increases with time as t 2 (if ho > 0), and the function 
e decreases as 1/t 2 , or vice- versa if ho < 0. In presence 
of the Lagrange multiplier, one can get exactly the same 
behavior as well as new solutions: e and V increase (for 



It is not difficult to express the Lagrange multiplier 
and the potential in terms of the other functions. In this 
way we get them by fixing the metric, the scalar field and 
the coefficient of the Gauss-Bonnet invariant. They are 



A = (46) 

= W-^-^-8* 2 !!-8;4(A-W), 



V = 



(47) 



2H + 3^ + £_ 8g2 ^_ 16g ^ + ff2) 



dt 2 



dt 



Using these expressions, it is easy to verify the validity 
of earlier solutions. 

Let us consider now the case when, as unknown func- 
tions do not appears the Lagrange multiplier and the 
function £((p). Combining the Friedmann Eqs. (1441) and 
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(flS")) . we obtain 



2 • 



= -zH + tp 2 + \(tp)tp 2 - 8H 



,,d 2 e{<p{t)) 



dt 2 



dt 



dt 

d ( H 2 de(ip(t)) \ 



^H + ip 2 + \(<p)<p 2 -8a- 
k at \ a at J 



Eq. (|48p can be solved with respect to e(<p(t)) as 



dt x 



8 J x ff(ti) 2 



(49) 



41 dU ( 2 



a(t 2 ) \k 



H{t 2 )+<p(t 2 ) 2 + \(<p(t 2 ))<p(t 2 y 



Combining Eqs. (|44l) and (|49|) . the scalar potential 



/t i± 
dk) x (50) 

x (^6 fa) + ^{hf + X{v(h))<p(h) 2 



Let us identify t with f((f) and iJ with g'(t) where 
/ and g are some unknown functions. The solution of 
the Friedmann equations is related to the existence and 
behavior of such functions. Then we consider a model 
where V(ip) and s(ip) can be expressed in terms of the 
two functions / and g as 



f'(, n Af>-9(f(<Pi)) 



3</(/M)e»<^» / dpi/'fai) 



AM 



-3 5 ' (f(<p))f*™» x 



(51) 



1 r , /Vje^OJ 

£ (W = o / X 



x / d(p 2 f'((p 2 )e 



1 r , /'(vi)es^ 1 )) 

' #1 - X 



5'(Vl) 2 



x / dtp 2 f'(tp 2 )e 



-g(f(v>2)) 



vi 



(M> ' (52) 

By choosing V(ip) and e(<p) as in Eqs. (fSTj) . one can easily 
find the following solution for Eqs. (|44l) and (|45[) that can 
be compared with results in [Hj|), that is 

= /-*(*) (t = m) , 

a = aoe»W (H = g'(t)) . (53) 

Hence any monotonic evolving cosmology, expressed 
by = <?'(</?) in the model @ with the potential (f5"Tj) . 
can be realized, including models exhibiting the transi- 
tion from non-phantom phase to phantom phase without 
introducing the phantom scalar field. However, we have 
to note that the approach could become problematic for 
H(t) non- monotonic as we will discuss in the next sub- 
section. 

As an example, let us consider the model (|28p. Hence, 



f{p) = t = he^°, /'(¥>) = — e 



(54) 



where prime is the time derivative (' = d/dt). Now we 
choose the Lagrange multiplier as in Eqs. (|4"2")h In this 
case, the potential and the function e iare composed by 
two terms: the first term is the value in the absence of a 
Lagrange multiplier. For e, it has the form 

t 2 (2h - <p 2 n 2 ) 2[p/lpo 
16fcg(l + h$)K 2 

If e 2Lp /vo = t 2 /tl, this expression is equal to (j3"5)) . The 
second term ([521 has the form 

2ip/vo (I6e 2 ^^ £0 ^(l + h )n 2 + (<p 2 K 2 - 2h )t 2 ) 
16h 2 (l + h ) K 2 

It also contains the first term with a different sign, and 
also the following expression 



(55) 
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which coincides with Eq. (l4"2"j) . This example is paradig- 
matic to show how the reconstruction scheme works for 
the Gauss-Bonnet gravity with Lagrange multiplier. 

A. The case of non-monotonic functions of time. 



where ip s = ipoto- We see that by choosing different types 
of functions A, we obtain different forms of V and e. 
It is not always easy to find the resulting integral: for 
example, choosing A = Ane^/^ , we obtain the solution 
in terms of the Euler function. However, it is easy to 
obtain solutions for the case A = An (ip/ipo) Xl ■ We have 



The above examples correspond to monotonic behav- 
iors in time for the Hubble parameter and the scalar field. 
However, from Eqs. (|5Tjl and (f52l) . one could run into sin- 
gularities emerging from non-monotonic behaviors (for 
example H — and <fi = for some values of t). If 
we consider the original Eqs. (|44[) and (|45|). we see that 
terms like 1/H and l/tp do not appear. The problem 
may occur when we recast V and e in terms of ip. To 
face the issue of non-monotonic behaviors and construct 
suitable examples, let us take into account the following 
form of the Hubble parameter 



H = ho 



1 



to-* 



(56) 



+ 



f 



32^2 

l5(p s K 2 



c 2 



2,„2 



(vo) 



3k 2 

Ai 



2tp 5 tp s 



25 4 

2 5 , V 

-ip (p s ln — 

5 ip s 



A 



(Ai - 1)A X (1 + Ai)(4 + Ai)(5 + Ai)(6 + Ai) 
x ( (<5 2 (A 1 -1)A 1 (4 + A 1 )(5 + A 1 )- 
- 2 Ws (A x - 1)(1 + Ai)(4 + A x )(6 + Ai) + 
+ ^A 1 (l + A 1 )(5 + A 1 )(6 + A 1 ))), 



(59) 



This means that the universe is in a non-phantom phase 
for t < to/ 2 and in a phantom phase for t > to/2. There is 
also a Big Rip singularity for t — to and a point where one 
can define an effective cosmological constant for t = to/2 
(£(t /2) = 0). 

One can take into account two different behaviors for 
the scalar field: 

1) ip = ipot, that is t = f(ip) = if/ifo and 

2) ip = (po In (t(to — t)), that is t = f(<p) — 

= \ (to ± y/tl - lev/vA . Clearly, expressing V and e 
by ip, uncertainties on the time evolution can come out 
while they are removed as soon as time-dependence in 
considered. 

Let us derive V and e for the first case. The general 
case contains special functions under the integrals and 
then we consider the specific case ho = 2. We have 



V 



4ip 2 (2ip + (p s ) 2 12<A) 2 W 

ip 2 (tp — <p s ) 2 K 2 

2ci(tp - ip s ) 2 + ipo 2 



, 2 ln^ 



2 ->- Ucpip s 



-f 



2{ip-ip s f 



(v°) 



An 



Qp - ^) 3 Ai (A 2 - 1) 
(-3(^ 2 ^(-l + Ai)A 1 (3 + A 1 )+ 
3 W 2 (-1 + A 1 )(1 + A 1 )(4+A 1 )- 
^A 1 (1 + A 1 )(5 + A 1 ) + (^ 3 A 1 (-1 



Uip s 



A?)) • (60) 



V = 



1 



32^ 2 
8ip 5 

15ip s K 2 

1 rf 



Clip° 



c 2 



Aip A ip a 2ip 2 ip 2 



2ip 5 ip s 
25 



4 



dip 



3k 2 3k 2 
//( yi )eg(/(yO) 

"7W ; 



-ip ip s In 

5 ip s 



■ip 2 (2ip + ips) 2 12( ^0 2 W 
ip 2 (p - ip s ) 2 K 2 

2c\{ip - ip s ) 2 + po 2 



, 2 ln^ 



(V 



- l^tpips 



X(ip) 



2{ v -p> s ) 2 
3g'(f(ip))e^» x 



\f'(<pi) 



(57) 



(58) 



This means that the delicate point is related to a suitable 
choice of A. 

Let us consider now the second case where the scalar 
field is a logarithmic function of time. We assume again 
ho = 2 and A = \oe~ 2(p ^ Vo . For this choice, we obtain 
solutions that do not contain special functions, that is 



ip 2 o /20ci< 7 



80i 2 V 7^ 2 
5 

24 : 



-ip 2 K 2 t 8 



2tH 



+ ipl^HHo + c 2 + 



5tH 2 



Aot^ 
2t\ 



3t 2 



Wo~W^ 

X t 5 Info - t] 
'o 



X Q t 5 ln[f] 
'o 



-A In [t - t] 



(61) 
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V 



<pI (At 3 + 4ot 2 t - mtl + 34) 
2t 2 (t - t ) 3 

At (-3ci« 2 < 3 + 2t {3t + t )) 

2n 2 t 2 (t-t ) 3 + 
p§A (-60t 5 + 30tH 



iot 3 ti - ibt 2 4 - 14 + 4) 



5t 4 {t- 

60ipl\ t(ln[t] - ln[* - *]) 



to) 4 t 3 



5(t 



to) 3 4 



In this case, we have written the expression as a func- 
tion of time, as we have the uncertainty noted above. 

In conclusion, it can be seen that the method works 
not only for monotonic functions, but also for some non- 
monotonic Hubble functions that have to be carefully 
considered. However problems may arise by converting 
the functions V and e in terms of the scalar field. In this 
process, uncertainty on the behavior can arise. The inter- 
est to derive suitable Hubble rate forms (not necessarily 
monotonic) is related to the necessity to obtain realistic 
dark energy models. In fact, H(t) decreasing with time 
is generally related to the evolution of matter dominated 
epochs and large scale structure formation. On the other 
hand, H (t) increasing with time can be related to infla- 
tionary and dark energy behaviors. Clearly, special care 
has to be devoted to singularities at finite which could 
come out, as above, in this reconstruction process. 



IV. DYNAMICAL SYSTEM ANALYSIS AND 
CRITICAL POINTS 

Let us rewrite the cosmological Friedmann equations 
in terms of an autonomous dynamical system. The ap- 
proach consists of two parts: the Friedmann Eqs. (JT]) and 
([5]) and equations specifying the model, for example (TT2")) 
and (|13l) . It is convenient to work with the new variables 

(63) 



H 



H 2 



The prime means the derivative with respect to In a. 

After some algebra, we obtain the following dynamical 
system 

= -4r + \x 2 + Xx 2 + z + 24y, 
K 2 

= (2u + 3 + -x 2 - z - 8y' - 8yu - 16y, 
K 2 



y' = y(b 2 x + -+2 f i), 



(64) 



z = b\zx — 2zn — 
A' = b 3 (X - X\)x , 



hAx 
-Jp-> 



where the first equation is a constraint for the other 4 
equations. For the sake of simplicity, we have considerd 
the case A = 0. As usual, the critical points at finite are 
given by the condition 



x = y = z = X = U . 
The various situations are 



(65) 



(62) 



• A: (x, y, z, X) = (0, 0, A) or = (0, 0, X x ) (two 
points - any A or A = Ai is constant). For this point 
/1 — and hence EoS parameter is equal to minus 
one (w = —1). Thus, it is a de Sitter space. In this 
point, the potential of the scalar field dominates. 

• R- (t „ 7 W-( . bl n 6« 2 (l+Ai)-b 1 (l+2A 1 ) \ \ 

• 0. [x,y,z,A)~( K 2 (1 4_ Al ) , u, 2^(1+17? ' M) 

b 2 

(A is constant). For this point [i = — 2^(1+17) anc ^ 



hence EoS parameter is equal: w 



3K a (l + Ai) 



1. 



Also in this point, the scalar field potential 
dominates. 



C: (x,y,z,X) = (±- 



s/6 



= , 0, 0, Ai) (A is con- 



",/k 2 (1+2Ai)' 

stant). For this point [i = — ^x^^ and hence EoS 

parameter is w = ixgx^ — 1- In this point the 
kinetic energy of the scalar field dominates. 

• D: For this point one gets complicated expressions 
for x and y (x ^ 0, y ^ 0), but z = and A = 
Ai = const Futhermore, one gets: a — — 

If A is not zero (we have to add one more equation 
H' = uH to the system (|M|). then the point A does not 
change. However, it appears a new point AI for which 
H = ±^S. This means that at least one more point 
with all of the variables not equal to zero exists. Follow- 
ing this analysis of critical points, it is straightforward to 
study the stability and attractor nature of cosmological 
solutions under discussion at finite and asymptotically. 

Finally we want to stress that the presence of the La- 
grange multiplier adds new critical points. Moreover, we 
can get a different value w for the same potential by the 
modification of the constraint. 



V. CONCLUSIONS 

We have examined a string-inspired effective theory 
of gravity containing Gauss-Bonnet invariant interact- 
ing with a scalar field where a Lagrange multiplier is 
inserted into the action. The Lagrange multiplier term 
can be related to the process of string compactification 
that leads to the effective 4D-action. Adopting a FRW 
metric, the corresponding cosmological Friedmann equa- 
tions ([7]) and ((HJ can be used to define the effective po- 
tential and the Lagrange multiplier (and the function e 
which is the coupling to the Gauss-Bonnet term). Choos- 
ing a form for the metric and the scalar field, one can 
easily approach the cosmological reconstruction in a sort 
of inverse scattering approach. In general, the presence 
of the Lagrange multiplier helps in the generation of 
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new cosmological solutions or in changing the features 
of some known solutions. This fact could be extremely 
relevant from an observational point of view because the 
Lagrange multipliers could be related to natural "pri- 
ors" in order to discriminate among concurrent models. 
Specifically, relating cosmographic observed parameters 
as {Hq, Qo, JOj so, f^Mj ^a} to some Lagrange multiplier 
could result a useful tool to restrict the range of models 
physically viable (see Ref . and references therein for a 
detailed discussion on this point). 

In the models that we have considered, the Lagrange 
multiplier can be interpreted as the presence of some dust 
fluid (matter without pressure) that affects the evolution 
of the cosmological system. In this case, the Lagrange 
multiplier can help to realize the transitions between 
matter dominance and dark energy era (and vice- versa) 
without imposing unnatural fine-tunings. The detailed 
study of the theory ^ with the exponential choice for 
the potential, the Lagrangian multiplier and the function 
e has been presented. Two forms of scalar field: tp ~ t 
and tp ~ In t are discussed. The comparison with the case 
of no Lagrange multiplier is considered. The presence of 
the Lagrange multiplier gives rise to new cosmological so- 
lutions due to the change in the dynamical system. For 
example, in the case of the scalar field with the logarithm 
form in time, the theory provides only one cosmologi- 
cal solution. In presence of the Lagrange multiplier, the 
number of cosmological solutions with different features 
increase. In particular, (ghost-free) phantom cosmology 
with canonical scalar easily emerges as solution. How- 
ever, the situation can be more complicated since the 
Hubble rate H is not necessarily a monotonic function of 
time. In some phantom models, dark energy grows with 
time and so does H at late times. On the other hand, H 
can decrease in matter-dominated eras so non-monotonic 
behaviors have to be considered. Above we gave also an 
example in this sense. 

The extension to other modified gravities (teleparallel, 
Horava-Lifshitz, F(G), non-local, etc Q) via the intro- 
duction of the Lagrange multiplier can be easily accom- 
plished in this framework. These topics will be discussed 
elsewhere. 



Appendix A: The Lagrange Multiplier Method and 
Noether Symmetries 

A more detailed discussion on the Lagrange Multiplier 
Method is necessary. A Lagrange multiplier is not intro- 
duced ad hoc in the dynamics but it is related to con- 
straints of the theory related to symmetries and conser- 
vation laws. Here we will sketch the Noether Symme- 
try Approach (see Ref.[18l| for a recent review) showing 
that the reduction of dynamics induced by the Lagrange 
multi plie rs is always related to the search for symme- 
tries |28l . [20( 1 . In particular, the form of the Lagrange 
multiplier derives from conservation laws and then it has 
always a physical meaning. 



We have used a FRW metric, so that the Einstein Eqs. 
© and © reduce to ordinary differential equations. The 
same equations can be derived from the action ((2j) where 
the Lagrangian becomes point-like after introducing in it 
a FRW metric. The configuration space is a minisuper- 
space where the Lagrangian coordinates are the scale fac- 
tor a and the scalar field ip, with the velocities a, ip [Isj . 
In this case, the dimension of the space is 2 but there 
are cases in which the minisuperspace dimension can be 
larg er (for instance a Bianchi universe with anisotropies 
3C I or dynamical systems with more than one scalar field 

With these considerations in mind, let £(q l ,q l ) be a 
Lagrangian, independent of time and nondegenerate, i.e. 



dC 



det Hi j d = det 



d 2 C 



dq l dqi 



^0, (Al) 



where det Hij is the Hessian determinant. In standard 
problems of analytical mechanics C has the form 



£ = T(q,q)-y(q) 



(A2) 



where T and V are the kinetic and potential energy. T is 
a positive definite quadratic form in q. Associated with 
C is the energy function 



w - dC * 



L. 



(A3) 



which is the total energy T + V. It is worth noticing that 
C can be more complicated than (|A2j) and Ec is a con- 
stant of the motion called "energy" also in general cases. 
In the Lagrangian formalism, we have to consider only 
transformations which are point-transformations. Any 
invertible and smooth transformation of the positions 
Q l = Q l (q) induces a transformation of the velocities 
of the form 



(A4) 



The matrix J = ||9Q l /^'7"' ll ^ s the Jacobian of the 
transformation on the positions and it is assumed to be 
nonzero. The Jacobian J' of the induced transformation 
is easily derived and J ^ — > J' 7^ 0. A point transfor- 
mation Q l = Q l (q) can depend on one (or more than one) 
parameter. Let us assume that a point transformation 
depends on a parameter e, i.e. Q l — Q J (q, e), and that 
it gives rise to a one-parameter Lie group. For infinitesi- 
mal values of e, the transformation is then generated by a 
vector field: for instance, as well known, d/dx represents 
a translation along x axis, x(d/dy) — y(d/dx) is a rota- 
tion around z axis and so on. In general, an infinitesimal 
point transformation is represented by a generic vector 
field on Q 



X = a l (q) 



d_ 



(A5) 
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The induced transformation (|A4[) , considering also veloc- 
ities, is then represented by 

+ w- (A6) 

X c is called the "complete lift" of X [20]. From now 
on, we will drop the suffix c but clearly we refer to a 
complete lift. A function /(q, q) is invariant under a 
transformation represented by X if 

§=°. (AT) 

where Lx/ is the Lie derivative of /. In particular, if 
Lx£ = 0, X is said to be a symmetry for the dynamics 
derived by C. 

Let us consider now a Lagrangian £ and its Euler- 
Lagrange equations 



d dC _ dC 
dt dqi dq j 



= 0. 



(A8) 



Let us consider also a vector field of the form (|A6j) . Con- 
tracting (|A8|) with the a"s gives 



d d£ <9£ 

dt dq> dqi 



= 0. 



(A9) 



Being 

dt dq> dt 
from (IA9[) , we obtain 
d / 



d d£ d f i dC\ (doP\ dC , Ain , 

;i^ = ^( a ^jj-(^rj^j' (A10) 



(All) 



The immediate consequence is the Noether Theorem: If 
Ly^C = 0, i/ien i/ie function 



dq l 



(A12) 



is a constant of motion. 



Eq. (|A12[) can be expressed independently of coordinates 
as a contraction of X with Cartan one-form 



def <9£ , i 



dq l 



(A13) 



For a generic vector field Y = y'^/dx 1 , and one-form 
P = Pidx 1 , we have by definition i^P — y i Pi- Thus 
Eq. (|A12p can be written as 



iy.de = So 



(A14) 



Under a point-transformation, the vector field X be- 
comes 



X' = ( ix dQ fc )JL + (±U xd Q k )^ ' S 



dQ k \dt 



dQ k 



(A15) 



We see that X' is still the lift of a vector field defined on 
the "positions" space only. If X is a symmetry and we 
choose a point transformation in such a way that 



ixdQ 1 = 1 ; i*dQ l = i^l 



we get 



X' 
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dC 



dQ 1 ' 8Q 1 



. 



(A16) 



(A17) 



Thus Q 1 is a cyclic coordinate and the dynamics can 
be reduced according to well known procedures (20l . [2l| . 
The change of coordinates defined by (|A16[) is not unique. 
Usually a clever choice is very important. However, more 
than one symmetry can exist and the degree of dynamical 
reduction depends on the number of symmetries [l8j . 

Let us now assume that £ is a canonical Lagrangian 
(e.g. of the form (|A"2jO . As X is of the form (|A15j) . L x £ 
will be a homogeneous polynomial of second degree in 
the velocities plus a inhomogeneous term in the q % . Since 
such a polynomial has to be identically zero, each coeffi- 
cient must be independently zero. If n is the dimension of 
the configuration space, we get {1 + n(n + l)/2} partial 
differential equations (PDE). The system is overdeter- 
mined, therefore, if a solution exists, it will be expressed 
in terms of integration constants instead of boundary 
conditions. Usually, C is fixed and one can ask for so- 
lutions of a specific PDE system. In particular, a La- 
grangian containing some undefined functions (e.g. cou- 
plings cj (</?), e(ip) and potentials V(ip)) is a class of La- 
grangians where the single element is determined by the 
solution of the PDE systems. In other words, the Noether 
symmetry selects the functions which assign the model. 

The above discussion shows that, given a dynamical 
system, it is always possible to search for Noether sym- 
metries. If they exist, the dynamics can be specified and 
reduced by a change of variables since one (or more than 
one) coordinate becomes cyclic. If we are able to in- 
tegrate such a new dynamics, the problem is to invert 
the solution in order to get the evolution in the previous 
variables [28l |29| . Conversely, the existence of Noether 
symmetries is able to select the form of scalar field po- 
tentials and couplings. In other words, given a class of 
models, the symmetries are able to select some of them, 
which, in principle, are of physical interest. 

The Lagrange multiplier method is related to this ap- 
proach. In other words, imposing the existence of La- 
grange multipliers means to select models where sym- 
metries exist because the Lagrangian functions become 
canonical. The example below shows the strict relation 
between the Lagrange multipliers and the Noether sym- 
metries. 

Let us consider the simple case of f(R) gravity defined 
by the action 



S 



d 4 a 



(A18) 



where, as usual, R is the Ricci scalar. This case is very 
interesting since not only the coupling and the potential 
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are unspecified but the whole Lagrangian is not given 
a priori. Imposing, as above, the FRW metric, we can 
write 



S= C(a 7 a,R,R)dt , 



(A19) 



considering a and R as canonical variables. Such a posi- 
tion is arbitrary since R depends on a, d and a so (|A19[) 
is not a true canonical Lagrangian. To remove this am- 
biguity, one can use a Lagrange multiplier A and search 
for Noether symmetry related to it. We have 



S = 2tt 2 / dt x 

,3 



(A20) 



x U(R)a 3 -X 



R + 6 



2f(R)^ + f(R) 



da df3 
2a + a— +a— 
Oa OR 



(A29) 

+ a(3f"'(R) = 
= 0, (A30) 



3a [f(R) - Rf\R)\ - apRf\R) = , (A31) 



af(R)+af3f"(R) = 0. (A32) 



From (|A29[) . we have that a is a function of a only, if 
we want to avoid trivial cases (i.e. we want f"(R) ^ 0). 
The symmetry is given by the functions 



To determine A, we have to vary the action with respect 
to R. that is 



<*=— , £ = -2/3o4, f(R) = foR 3/ \ (A33) 



a 3^M SR _xSR = 0, 



dR 



(A21) 



from which we get, like in Sec. II, the functional form of 
the Lagrange multiplier, that is 



A = a 3 f'(R) 



(A22) 



Substituting into (|A20j) and integrating by parts, we ob- 
tain the point-like Lagrangian 

C = a 3 [f{R)- Rf(R)]+6a 2 af'(R) + 

+ Qa 2 aRf"{R)~akf{R). (A23) 

Then the cosmological equations of motion are 



my 



2{-)f(R)"R + 



(A24) 



+ f{R)"R + f{R)"'R 2 - - [Rf(R)> + f(R)\ = , 



and 



which solve the above system; /3o and /o are integration 
constants. The new induced variables can be 



w = a R, z = — , 
from which the Lagrangian (|A23|) becomes 



which can be rewritten in the form 



£ = 9p zy - 9fcy - -y 3 , 



(A34) 



~ Q 7,111 \ i 

C= ^ ^-9kV^--V^, (A35) 
2 %Jw 2 



(A36) 



using y — y/w. Dynamics is then described from the 
equations 



y = , from which y = y = S : 



(A37) 



, a a k 
i?=-6(- + — + — 

a a z a z 



(A25) 



9/3 z + 9k + -y 2 = ; 



(A38) 



where, as above, the Lagrange multiplier gives one of the 
equations of motion. The energy constraint is 



6d 2 af'(R) - a 3 [f(R) - Rf'(R)} 
6a 2 dRf"(R) + akf'(R) = 0. 



(A26) 



The symmetry generator is defined on the tangent bundle 
TQ{a 1 d, R, R) and it is 



X = aM A +/3 (a,*)A 



da d df3 d 
~dtd& + dtdR' { ' 



while the Noether condition L^C — produces the sys- 
tem 



„ da 
a + 2a— 
oa 



af"(R) 



oa 



= 0. 



(A28) 



9/3 yi + 9ky + —y = , 



whose solutions are 



v(t) = yot + yo ■ 



z (t) = c A r + c 3 t 3 + c 2 t 2 + at + c . 



ith 



C4 



C2 



Vo _ * 
12/? 2 • 



yoyo 
3/3 



(A39) 



(A40) 



(A41) 



Cl = z , c = Zq . 



(A42) 
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The energy condition (|A39I) gives the relation among the 
initial data [29] . Going back to the physical variables, we 
have 



a(t) = ±^d 4 t 4 + d 3 t 3 + d 2 t 2 + dit + d , (A43) 

where the dj's are the c,'s multiplied by 2/?o- The R 
variable is actually the Lagrange multiplier which gives 



R 



(yot + yof 



d 4 t 4 + d 3 t 3 + d 2 t 2 + dit + d 



(A44) 



and then the cosmological equations (|A24|) . (|A25|) and 
(|A26[> arc fully satisfied. In conclusion, by imposing the 
Lagrange multiplier we get a canonical dynamics and the 
consequent existence of the Noether symmetry allows its 



integration. In summary, the dynamical system results 
both canonical and integrable. More complicated cases 
are discussed in [H, [l^| . 
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